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We investigate collective effects of an ensemble of biased two level systems interact- 
ing with a bosonic bath in the strong coupling regime. The two level systems are 
described by a large pseudo-spin J. An equation for the expectation value M(t) 
of the ^-component of the pseudo spin is derived and solved numerically for an 
ohmic bath at T = 0. In case of a large cut-off frequency of the spectral function, a 
Markov approximation is justified and an analytical solution is presented. We find 
that M(t) relaxes towards a highly correlated state with maximum value ±J for 
large times. However, this relaxation is extremely slow for most parameter values 
so as if the system was "frozen in" by interaction with the bosonic bath. 

1 Introduction 

The spin-boson model deals with a biased twoJevel system, described as a spin |, 
under the influence of a bosonic environment .□ It is applied to different physical 
systems such as, e. g., two level systems in amorphous materials, or double quantum 
dots interacting with phononso'Em or photons. Further examples include a magnetic 
flux in a SQUID which originally motivated the spin-boson model.El 

In this contribution, we study the more general case of many identical two level 
systems (or larger spins) that couple to a bosonic field. Examples are glasses or 
arrays of coupled dots. We focus on collective effects and hence, investigate the 
dynamics of the total spin with projection J z — Y] f erf. This pseudo spin gives 
information about the degree of polarization, i. e. the number of two level systems 
to be found in each of the two states. We concentrate on the strong coupling regime, 
when perturbation theory is not sufficient to take into account the interaction with 
the environment. 

2 The large spin boson model 

The large spin-boson system is described by the Hamiltonian 

H — eJ z + 2T C J X + J z ^2 7 g (a_ 9 + a+) + ^ u q a+a 9 , (1) 

9 9 

which is formally equivalent to the spin-boson Hamiltonian. The only difference 
is that J, obeying the spin algebra [Ji, Jj] = SijkiJk, is not restricted to a spin | 
but can take any value, e and T c are the energy difference and the tunnel rate, 
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respectively, between the two states of the two level systems. The environment is 
given by the bosonic field operators a q , eigenenergies uj q (h — 1) and interaction 
coefficients 7 9 for the mode q. To take into account arbitrarily strong couplings 
between the large spin and the bosonic bath, a polaron transformation is applied 
to the Hamiltonian, 

H = e aJ 'He~' yJ ' = e J z - a J z + u q a+a q + T c ( J+ X + J_ X + ), (2) 

<? 

1 2 

with a :=V ^-(-a- g + at), X := e CT , a: =YL^L, ( 3 ) 

As a result of the polaron transformation the interaction term between the spin and 
the environment disappears whereas the tunnel term 2T C J X becomes dependent on 
the bosonic field. Furthermore, a new term —aJ z appears in the transformed 
Hamiltonian which will later turn out to be of importance for the long time be- 
haviour of the large spin. For a single two level system, i.e. J = i this term does 
not play any role since J 2 is constant in that case. 

Next, the transformed Hamiltonian (|J) is used to derive a master equation for 
the density matrix p(t) of the spin. In second order Born approximation,!?! its time 
evolution is given by 

'm = -J dt 1 Tr Res {[H T (t), [H T (t'),p(t') ® Ra]]} , (4) 

where i/r(i) is the the tunnel part of the Hamiltonian in interaction picture. In 
this equation, only a product of the spin density matrix p and the density matrix 
of the bosonic bath Rq (which is assumed to be in thermal equilibrium at all times) 
appears. This makes it possible to trace over the reservoir's degrees of freedom 
without knowing the full density matrix of the combined system of spin and envi- 
ronment at all times. Note that pit') is evaluated at times t' in the integral in Eq. 

, thus describing a non-Markovian process. Although we used the second order 
Born approximation, the interaction between spin and bath is taken into consid- 
eration to arbitrary high orders due to the polaron transformation. A closed set 
of equations for the diagonal elements pm,m in the basis | JM) of eigenstates of 
J 2 and J z is derived from the master equation (|J) under the assumption that the 
elements pm,m±2 can be neglected. These correspond to simultaneous spin flips of 
two spins and hence two boson processes. In the semiclassical limit of many two 
level systems and therewith large spin J,Q the variance of J z can be neglected and 
we find for the expectation value M(t) = (J z ) t = Em ^Pm,m(^)> 

M(t) = -2T C 2 fdt'(j+M(t'))(j-M(t') + i) Re|e i(e - a(2M( * ,)-1))( "' ) C(t-t / )} 

- (J+M(t') + I)(j-M(t')) Re{e- 4(£ - Q ( 2M (*') +1 »("')C7(t-t')}- 

(5) 
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3 The bosonic bath 



The correlation function 

C(t-t') = Trues {Ro X t X+} (6) 

describes the influence of the environment on the large spin in Eq. (^). It acts as 
a memory that returns information about the spin at time t' back to the spin at 
time t after the delay t—t'. As usual, we assume that the spectral function J (to), 

JH = ^| 7? | 2 <5( W -^), (7) 

exhibits a power law behaviour J(uj) = g o; s e _ "/" c with an exponential cut-off 
at frequency w c , and g the effective interaction strength. Then, the correlation 
function can be calculated analytically. In the case of an ohmic bath (s = 1) and 
at temperature T — 0, the correlation function reads 

C(t) = (l + iuj c ty 9 . (8) 

For the limit of an infinite cut-off frequency lo c , we find 

Jdm^c Re{e^C(t)} = 0(^) S(t). (9) 

Thus, the Markov approximation is justified for large cut-off frequencies where 
analytic solution of Eq. (|5|) becomes possible. For g = 1 the inverse function 
t = t(M) is 

t(M) = { ((2J+1)rrl (/(M) " /(Mo)) ' Mo < -5 + m 

with T = e-e/^-^nTZ/wc, t = e^^-^nT? /uj c and 

f(x) = er 2{J+1) Ei(-2x+2(J+l)) - e 2J Ei(-2x — 2 J). (11) 

Ei(x) is the exponential-integral function and a = guj c . Initial values Mq in the 
interval Mq G [—1/2 + e/(2w c ),l/2 + e/(2u! c )] are not considered here, as this 
interval is of width 1 and hence contains only one possible value of Mq. 

4 Results 

We solved the integro-differential equation (||) for the expectation value M(t) nu- 
merically. This enables us to check the Markov approximation which was applied 
to derive the analytical solution ([l0|). In case the assumptions for the derivation 
of Eq. ( p^| ) are not fulfilled (g ^ 1 or small w c ), results can only be obtained by 
the numerical solution. Fig. [l].a illustrates that the inversion M(t) approaches + J 
or — J for large times depending on the initial value of Mq. If Mq > e/2ui c , we 
find that for large times t, M(t) relaxes towards a state with macroscopic complete 
inversion + J which means that all spins are up in spite of e > 0. This behaviour is 
due to the spin part eJ z — aJ 2 of the transformed Hamiltonian (0), the eigenvalues 
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Figure 1. a shows M(t) for eight two level systems, i. e, J = 4. with e = T c , u c = 20T c and g = l 
for initial values Mo=±l. The analytical solution (EqJlQ|) and the numerical solution of Eq.(^) 
coincide, indicating that the Markov approximation is justified. The inset shows the spin system 
eigenenergies as explained in the text, b: Analytical solution for M(t) for J = 10, e = 440T C , 
lu c = 20T c and g = l. After an initial phase, the decay is strongly damped as can be seen from the 
inset, which shows M(t) for larger times. 



of which lie on an inverted parabola with two minima at the extreme points ±J 
(inset of Fig. |l|a). Therefore, for very strong coupling to the bosonic field, the total 
system minimizes its energy by relaxing towards either to the global (M = — J) 
or local (M = + J) minimum. This situation is analogous to the emergence of a 
macroscopically polarized state in the strong coupling limit of the one-mode su- 
perradiance model.Lru However, the dynamics of this relaxation is extremely slow 
if M{t) is not close to e/2uj c as can be seen from Fig. [|.b (the parameters were 
chosen such that M{t) decays to — J even for the initial value Mo = +</)• Thus, the 
system seems to be "frozen in" by the interaction with the bosonic bath for most 
parameter values. 
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